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Abstract. We have incorporated time-dependent terms into the statistical and radiative equilibrium 
calculations of the non-LTE line-blanketed radiative transfer code CMFGEN. To illustrate the 
significant improvements in spectral fitting achieved for photospheric phase Type II SN, and to 
document the effects associated with time dependence, we model the outer 6.1 M of ejecta of a 
BSG/RSG progenitor star. Hopping by 3-day increments, we compute the UV to near-IR spectral 
evolution for both continuum and lines, from the fully ionized conditions at one week to the partially 
recombined conditions at 6 weeks after the explosion. We confirm the importance of allowing 
for time-dependence in the modeling of Type-II SN, as recently discussed by Utrobin & Chugai 
for SN1987A. However unlike Utrobin & Chugai, who treated the radiation field in a core-halo 
approximation and assumed the Sobolev approximation for line formation, we allow for the full 
interaction between the radiation field and level populations, and study the effects on the full 
spectrum. At the hydrogen-recombination epoch, Hi lines and NaD are considerably stronger and 
broader than in equivalent steady-state models, while Call is weakened. Former successes of steady- 
state CMFGEN models are unaffected, while former discrepancies are cured. Time dependence 
affects all lines, while the continuum, from the UV to the optical, changes only moderately. We 
identify two key effects: First, time dependence together with the energy gain through changes in 
ionization and excitation lead to an over-ionization in the vicinity of the photosphere, dramatically 
affecting line optical depths and profiles. Second, the ionization is frozen-in at large radii/velocities. 
This stems solely from the time-scale contrast between recombination and expansion and will occur, 
modulo non-thermal excitation effects, in all SN types. The importance of this effect on spectral 
analyses, across SN types and epochs, remains to be determined. 



INTRODUCTION 

Because of radiative cooling and the fast expansion of the exploding mantle of the 
progenitor massive star, photospheric-phase Type II supernova (SN) spectra evolve 
rapidly. At shock breakout the spectral energy distribution (SED) peaks in the far-UV 
and X-rays. Subsequently the SED centers in the UV, then in the optical, and later 
in the infrared until the object disappears altogether into oblivion after a few years. 
From a radiative transfer perspective, the cooling induces a recombination to lower 
ionization stages that impose a strong blanketing effect on the energy distribution. A 
few days after explosion, the UV/optical/IR spectrum shows a UV-dominant SED, with 
Balmer/Paschen lines of Hydrogen, Hel5875A, Hei 10830A, Nn4600A, Nll5400A, 
and a few isolated resonance lines such as Mgll 2802A and Aim 1859A. About a week 
later, hydrogen recombination kicks in, tied to a sudden enhancement in line-blanketing 
due to the switch from Fein to Fell. Dessart and Hillier yj] successfully modeled 
these epochs in the evolution of SN1999em, using steady-state non-LTE CMFGEN 



models, solar metallicity, and a power-law density decline with an exponent of ten. In 
particular, Hel lines at early times were fitted without invoking any additional non- 
thermal excitation due to 56 Ni, as required by Baron et al. M\ m their analysis of 
SN1993W. In Dessart and Hillier yl, the spectroscopic analysis deliberately focused 
on the first 40 days since at such an epoch, we encountered growing difficulties in 
reproducing Ha, H/3, and the near-IR Paschen lines [3]. While we could reproduce 
the overall SED very well, line profiles were sizably narrower than observed at late 
times, suggesting that line formation regions predicted by CMFGEN were confined to 
velocities/radii that were too small. CMFGEN failed to reproduce observations, despite 
all attempts, whenever hydrogen started recombining at and above the photosphere. 

The Ha problem has not been clearly emphasized in the literature, where one can 
see a great disparity in model atmosphere assumptions and agreement between theoreti- 
cal predictions and observations, but no direct links to physical/numerical issues. In the 
eighties and early nineties, the recognized importance of non-LTE effects confronted 
the strong limitations of computer technology, so that only a few species were treated in 
non-LTE (i.e., usually hydrogen and helium), while the metals responsible for line blan- 
keting were treated in LTE. Eastman and Kirshner [4] followed this approach to model 
the first ten days of SN1987A, and thus did not encounter obvious difficulties with hy- 
drogen lines. Schmutz et al. o], using an approximate non-LTE technique, had, on the 
contrary, great difficulty reproducing any of the Balmer lines, suggesting clumping as 
the culprit. Hoflich reproduced the SN1987A spectral evolution and the hydrogen 
lines, over many months, using a large "turbulent" velocity. In the more sophisticated 
non-LTE CMFGEN flU models presented in Dessart and Hillier 1 1, DH06], we found 
that decreasing the turbulent velocity weakens line-blanketing effects and increased, al- 
though only modestly, the strength of hydrogen lines. However, beyond 40 days after 
explosion, this tuning had no longer any important influence on the hydrogen lines. The 
non-LTE model atmosphere code PHOENIX reproduces the Balmer lines well |2|,|9j], but 
this may stem from their adoption of non-thermal ionization/excitation due to 56 Ni at the 
photosphere just a few days after explosion, in mass shells moving at >10000km s _1 in 
SN1987A |9] or ~9000km s _1 m SN1993W H. Hydrodynamical simulations of core- 
collapse SN predict that 56 Ni has velocities of at most ~4000km s _1 , the nickel fin- 
gers being strongly decelerated at the H/He interface ifiol [Till . The magnitude of this 
disagreement extends far beyond the uncertainties of explosion models and suggests a 
genuine incompatibility. 

Interestingly, Hi lines remain strong for months in all Type II SN, e.g. the "peculiar" 
SN1987A, the "plateau" SN1999em, and the "low-luminosity" SN1999br. CMFGEN 
models computed for these objects were unable to reproduce Balmer lines after ~4 days 
in SN1987A, -40 days in SN1999em, and ~20 days in SN1999br, all coincident with 
hydrogen recombination in the outflow. These three SN boast very different envelope 
structures, ejecta velocities, and light curves. SN1999br even synthesized an order 
of magnitude less 56 Ni than average LI 2]. The only common property between these 
objects, which is connected to the Ha problem, is the recombination of the ejecta to a 
lower ionization state at the corresponding epoch. 

Recently, Utrobin and Chugai lfl3l UC05] proposed that the effect of time- 
dependence, and the energy associated with changes in ionization/excitation, lead 
to a strong Ha line profile in SN1987A during the recombination epoch. They also 



found that barium lines were affected, and that, with their more consistent approach, 
Ball 6142A could be fitted using the LMC metallicity value. The steady-state models of 
Mazzali et al. IU4I1 supported instead an abundance enhancement of five. Time depen- 
dence has been invoked in the past by Fransson and Kozma 11 511 to explain the late-time 
light curve of SN1987A, and the theoretical study of Pinto and Eastman lfl6l [l7ll 
showed that time dependence in the radiation field had a critical impact on the radiative 
transfer in Type las. Pinto and Eastman IU6LI17I1 . however, treat the material in LTE, in 
other words, they neglect any non-LTE and time-dependent terms in the rate equations 
(i.e. they do not solve the rate equations), focusing instead on the time-dependent 
diffusion of photons through an optically thick Type la SN ejecta. Similarly, Kasen 
et al. O neglect such explicit time dependence in the rate equations, computing the 
ionization and excitation state of the medium in LTE. Thus, although there is at present 
growing interest in accounting for time dependence in the radiation field, the often-used 
expedient of LTE, to maintain CPU costs low, has forced the neglect of both non-LTE 
and time-dependence in the level populations. One exception to these time-dependent 
LTE approaches is the work of Hoflich lfl9ll . who treats time-dependence in the rate 
equations but, to our knowledge, has not discussed the associated effects on Type II SN 
spectra. In the present and forthcoming studies, we investigate thoroughly the effects 
of time-dependence in the rate equations and their impact on inferred ejecta properties. 
Time dependence in the radiation field is accounted for by ensuring that the emergent 
light from the model matches the observed flux from the Type II SN, using SN1999em 
as our reference SN yj]. 

Here, we report the salient features of a time-dependent non-LTE CMFGEN simula- 
tion covering the evolution of a Type II SN from one to six weeks after core collapse. 
We confirm the results of UC05 that time-dependence induces an over-ionization of the 
recombining ejecta and that it solves the Ha problem. However unlike UC05, we self- 
consistently solve the radiation transfer equation — the coupling between the level pop- 
ulations and the radiation field is calculated and fully allowed for. Our study of the full 
spectrum predicts that all lines are significantly affected. The ionization of the ejecta and 
its evolution are so strongly modified that we pred ict even Hel 10830A many weeks after 
explosion. This feature has been observed 0. 12011 . but never understood, and epitomizes 
the difficulty of line identification, abundance determinations, and chemical stratifica- 
tion in SN ejecta. In the next section, we present our baseline model, before moving 
onto the presentation and discussion of our results. 



MODEL 

CMFGEN 0, li solves the radiative transfer equation in the comoving frame, subject 
to the constraints of radiative and statistical equilibrium. It assumes spherical symmetry. 
No explicit time-dependence in the radiation field is considered, i.e. light propagates 
instantaneously and thus is not subject to travel-time delays [21], and we still ignore 
relativistic effects. Both affect the energy transport (luminosity) at depth lfl6l[l7ll . How- 
ever, unlike the version of CMFGEN used in il|22l[!J], we now include explicitly time 
dependence in the rate equations that are solved for level populations, as well as in the 



energy equation by accounting for changes in ionization and excitation. Our treatment, 
which follows directly from Mihalas and Weibel Mihalas II23I1 . will be presented in detail 
in a follow-up paper (Hillier & Dessart 2006, in prep.). 

We adopt a different procedure to model SN spectra from that employed previously. 
Rather than modeling, at each time, the region in the vicinity of the photosphere (i.e., 
from TR 0Sse i an d of 30), we model the time evolution of a larger portion of the ejecta. 
As a starting point for our baseline model, we mapped onto the CMFGEN grid an 
envelope with the following properties: Rq — 2.524 x 10 14 cm (Rq is the base radius), 
vo = 3470km s _1 (vo is the base velocity), R^po = 6.781 x 10 33 g (po is the base density), 
p(r) = po(Ro/r) 10 (r is the radius), and Lq = 3.6 x lO u L (Lo is the base luminosity). 
The maximum radius is 16 times the value of the base radius Rq. These initial conditions 
correspond to a base Rosseland optical depth of 4000, a total mass of 6.1 M Q , and, 
assuming homology, a time after explosion of 8.4 days. We adopt a fixed composition, 
i.e. chemical homogeneity, with values consistent with the red-sugergiant progenitor 
of SN1999em yJ,|24J]. Adopted abundances, given by number, are: H/He = 5, C/He = 
0.00017, N/He = 0.0068, O/He = 0.0001, and all metal species are at the solar value. 
Compared to the model atom used in DH06, we omit Cl, neon, aluminium, sulfur, 
chromium, manganese, and cobalt. 

An initial steady-state model is then computed with the above parameters and it gives 
a spectrum that fits the observations of SN1999em at 7 days after explosion [81 Fig. 1], 
with a fully-ionized ejecta, little blanketing in the UV, and a dominance of hydrogen 
Balmer lines and Hel5875A in the optical range. We then hop in time by steps of 
At = 3 days, i.e., all mass shells m in the initial model are evolved from their position 
r(m) to a new radius r'(m), with r'(m) = r(m) +vo(m)At, a new density p' = p(r/r') 3 , 
but the same velocity. In practice, for each new model, CMFGEN adapts its 80-point 
grid to cover each optical-depth decade with at least five points. 

The luminosity, set at the base, is adjusted at each epoch so that the emergent lumi- 
nosity matches the inferred values for SN1999em (DH06), i.e. from ~10 9 L one week 
to ~ 10 8 L 50 days after explosion. Despite the fixed composition at all times and the 
neglect of Cl, both of which affecting the late time appearance of Type II spectra (DH06) 
we support our thesis by showing in the next section a representative match between the 
spectrum of our baseline model and the optical observations of SN1999em. Detailed 
analyzes will be presented in a future paper. 



RESULTS AND DISCUSSION 

To highlight the salient features of time-dependent CMFGEN models, we select the last 
baseline model computed in the sequence, corresponding to 48.7 days after explosion, 
and compare its properties to the same model, but computed by assuming steady-state. 
Fig.[T]shows the radial variation of the ionization fraction of hydrogen and helium for the 
time-dependent and steady-state models (left and middle panels). Time dependence sys- 
tematically reduces the effective recombination rates, with the greatest effect occurring 
away from the photosphere where the densities are lowest and the velocities highest. In 
those regions and at 48.7 days, the ejecta preserves the high ionization it had at the start 



of the time sequence, at ~8 days after explosion, exhibiting the prevalence of ionized 
hydrogen and helium. Many other species are also affected, often showing comparable 
ionization profiles (i.e., a reduction of their degree of ionization near the photosphere 
and a growth of their degree of ionization at larger distances). Further, because of the 
enhanced H ionization, the electron density does not fall as steeply in the time-dependent 
CMFGEN model (Fig.[T]). The electron temperature is strongly affected in the optically- 
thick layers, which remain hotter due to the reduced radiative cooling efficiency in the 
time-dependent case. 

Thus, the effect of time dependence modifies the optical-depth of both continua and 
lines. The former occurs since more free electrons are available to scatter continuum 
photons. Further, certain ions are "over- abundant" by several orders of magnitude com- 
pared to the steady-state model. In the top panel of Fig. 01 although not tailored specifi- 
cally to match an observation, the reddened (E(B — V) = 0. 1) spectrum of our baseline 
model (black) is overlaid on the optical spectrum of SN1999em observed on the 5th of 
Dec 1999 JH, blue], i.e. ~45 days after explosion (DH06). Note the strong Ha line in 
this time-dependent CMFGEN synthetic spectrum, in aggreement with the observation. 
In the lower panels, we present synthetic UV, optical, and near-IR spectra for the time- 
dependent and steady-state (red) CMFGEN models, normalized at 12000A. To highlight 
the effect on individual species, we also compare synthetic spectra computed by allow- 
ing all continuum processes, but accounting only for bound-bound transitions of Fell 
and Hi. Note the large effect on Balmer and Paschen lines, which appear stronger in 
both absorption and emission and are broader. Fell lines are stronger in the steady state 
model — stronger in absorption in the UV and stronger in emission in the optical. Sim- 
ilarly, in the lower panel, we zoom in on specific synthetic lines, accounting only for 
bound-bound transitions of Nal, Call, or Hel. 

In the time-dependent case, Nal 5895A strengthens and broadens. This occurs because 
Na + is the dominant ionization state in both models, but the increased electron density 
in the time-dependent model leads to an increase in the neutral Na fraction. Conversely, 
Call 8500A weakens, primarily because the relative fractions of Ca + and Ca ++ change 
significantly. Hel 1.083 /im is present only in the time-dependent case, as a broad flat- 
topped profile (overlapping with Py in the full spectrum; top panel of Fig.O, with a blue 
absorption maximum at 1.038 /im (equivalent to —13000km s _1 ). 

The implications of time dependence are numerous and far reaching. We find that 
the effects of time dependence on level populations and on the energy budget allow the 
best of two worlds. Previous successes of CMFGEN at fitting the overall spectral energy 
distribution of photospheric phase Type II SN are left unchanged. The outstanding prob- 
lem of the Ha line strength and the overall line profile widths are however solved. We, 
thus, confirm the importance of time-dependence effects, as suggested by UC05, even at 
early times. This provides a more physical and less dramatic explanation of the appar- 
ent ionization/excitation seen in photospheric-phase Type II SN spectra than invoking 
non-thermal ionization/excitation by 56 Ni/ 56 Co in ejecta regions moving sometimes at 
~10000kms -1 . 

The change in ionization is large. Time dependence affects line profiles in a way 
that can mimic changes in chemical abundance in steady-state models, with dire con- 
sequences on the quantitative assessment of the chemical composition of SN outflows. 
Even more extreme is the case of Hel, absent in our steady state models after ~2 weeks, 




FIGURE 1. Left: Radial variation of the ionization fraction of Hydrogen (Hli: solid line; Hi: Broken 
line) for the models with time-dependence (black) and without (red), and at 48.7 days after explosion, 
but otherwise having identical parameters. Note the overionization of hydrogen at large radii, not present 
in its neutral state despite the very late time, akin to a frozen-in ionization. Symbols give the positions 
of the CMFGEN adaptive-grid points, for each model but for only one curve. Note the presence of the 
well-resolved recombination front near the base, deeper in for the steady-state CMFGEN model. Middle: 
Same as left, but for the ionization fraction of helium (Hell: solid line; Hel: Broken line). Right: Same as 
left, but for electron density and temperature. 
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FIGURE 2. Comparison of the synthetic spectrum for the time-dependent CMFGEN model (black 
curve) at 48.7 days after explosion with Keck observations of SN1999em II25II on Dec. 5th 1999 (turquoise 
curve), at ~45 days after explosion (DH06). The model has been reddened with E(B — V) =0.1. 



but present in equivalent time-dependent models even after 6 weeks, without any en- 
hancement in abundance from the initial red-supergiant composition. An unidentified 
feature at 1.03jUm was observed by Hamuy et al. yD in SN1999em at 4 weeks after 
explosion but not reproduced by DH06 with steady-state CMFGEN models; with time 
dependence effects, it can now be identified as Hel 1.083 /im. It also seems to be present 
in SN1987A a few weeks after the explosion lEoh . It has been inferred in Type Ib/c SN 
lEfill . produced by non-thermal excitation ll27ll . and even suggested in Type la SN lEsIl , 
although Mgll is expected to contribute, perhaps dominantly, to the feature 112911 . In the 
future, we will investigate to what extent time dependence can be distinguished from 
non-thermal excitation, and if, together with the accurate treatment of non-LTE effects, 
they can help resolve the issue of the presence of helium in SN spectra. 

Our investigation suggests that to model the emergent light from a SN at a given 




FIGURE 3. Top: Comparison of the synthetic spectra for the time-dependent CMFGEN model (black 
curve) and the equivalent steady-state CMFGEN model (red curve), both at 48.7 days after explosion. 
Accounting for time dependence leads to a comparatively higher ionization in the ejecta, sustained line 
optical depth at large radii/velocities, epitomized by Ha. Middle: Rectified spectra (shifted vertically) 
obtained by accounting solely for bound-bound transitions of Fell (lower curves; line spectrum magnified 
by a factor of four beyond 5500A for better rendering) and Hi (upper curves). Bottom: Same as above, 
but this time zooming-in on specific spectral regions, i.e., Nal5895A (left), Call8500A (middle), and 
Heil.083/im (right). 



time, one needs to have the knowledge of the physical conditions in the SN envelope 
at previous times, all the way to the epoch when the ejecta was fully ionized. This 
calls for a change of strategy in SN spectroscopic modeling, and casts doubt on results 
using approaches that ignore the history of the SN ejecta. While the energy gain through 
hydrogen recombination in the outflow is a feature of Type II SN ejecta only, the frozen- 
in ionization at large distances/velocities stems from the order-of-magnitude larger 
recombination time-scale compared to expansion time-scale. This fact is not influenced 
by the energetics of recombination — it merely demonstrates that SN ejecta are so 



tenuous and so fast that the ionization balance can be considerably affected by time 
dependence. This phenomenon should be a ubiquitous feature of SN ejecta and, at the 
very least, deserves considerably more attention that it has so far received. 

Accounting for time dependence in CMFGEN opens new horizons. In particular, 
it is now sensible to model Type II SN throughout the photospheric phase. We will 
incorporate in CMFGEN the energy deposition due to unstable isotopes, to cover the 
nebular phase, and investigate a wider range of SN explosions, from core-collapse to 
thermonuclear runaways. We will also include relativistic effects, and time dependence 
for the radiation transport, to improve the energy transport at depth. We will finally 
assess the effects of time dependence on the correction factors used in the Expanding 
Photosphere Method lEil. 
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